Bose-Einstein condensation in arbitrarily shaped cavities 
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We discuss the phenomenon of Bose-Einstein condensation 
of an ideal non-relativistic Bose gas in an arbitrarily shaped 
cavity. The influence of the finite extension of the cavity on 
all thermodynamical quantities, especially on the critical tem- 
perature of the system, is considered. We use two main meth- 
ods which are shown to be equivalent. The first deals with 
the partition function as a sum over energy levels and uses 
a Mellin-Barnes integral representation to extract an asymp- 
totic formula. The second method converts the sum over the 
energy levels to an integral with a suitable density of states 
factor obtained from spectral analysis. The application to 
some simple cavities is discussed. 



I. INTRODUCTION 

It is now well over seventy years since the phenomenon 
referred to as Bose-Einstein condensation (BEC) was first 
predicted jj]]2). Until recently the best experimental ev- 
idence that BEC could occur in a real physical system 
was liquid helium, as suggested originally by London j| . 
More recently it was suggested [[||| that BEC could oc- 
cur for excitons in certain types of non-metallic crystals 
(such as CuCl for example). There is now good evidence 
for this in a number of experiments |l. However, the 
most exciting experimental evidence for BEC has come 
from the observations of very cold alkali gases. BEC 
has now been observed to occur in gases of rubidium |Q , 
lithium H and sodium 0], In these experiments BEC 
has been realized in magnetic traps and laser traps, in 
small volumes far away from the infinite volume limit. 
Although nowadays a few million particles remain in the 
trap, in these early experiments only up to 2000 particles 
remained there so that it was by no means clear if calcula- 
tions done in the thermodynamical limit are an accurate 
treatment when considering this experimental situation. 
(For a more detailed analysis of this see Refs. |L0] [L5|.) 
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This motivates us to consider the related problem of the 
quantum statistics of a finite number of particles in an 
arbitrarily shaped cavity. 

When thinking about quantum statistics in finite vol- 
umes there appears an immediate difference to the ther- 
modynamic limit. Quantum particles confined to a finite 
volume of arbitrary size inevitably have a nonzero kinetic 
energy, hence they must exert some pressure. Taking 
/^m-sized cavities and the parameters typical for liquid 
helium it can be seen that the zero-point pressure is of 
measurable magnitude E3j. But not only the pressure 
but also all other thermodynamical quantities may devi- 
ate substantially from the results obtained in the thermo- 
dynamical limit once the volume and the particle number 
get small enough. In this article we want to give a sys- 
tematic treatment of this problem. 

In doing this we are going to extend a recent careful 
and beautiful study of BEC in a cubic volume where the 
gas was supposed to fulfill Dirichlet boundary conditions 
at the boundary of the cavity Jll|. The authors there 
employed an approach where the sum over the discrete 
spectrum was replaced by an integral with an appropri- 
ate density of states. This density was obtained by the 
knowledge of the energy spectrum in this example. 

In this article we are going to treat an arbitrarily 
shaped cavity. The methods we are going to employ are 
the so called heat-kernel techniques extensively used in 
finite temperature relativistic quantum field theory start- 
ing with the work of Dowker, Kennedy fl7| . However, in 
these considerations the stress was more on the influence 
of gravitational fields on the quantum statistics (see also 
p8[ ) . Recently also the influence of boundaries was con- 
sidered in the context of the Casimir effect |19j2^| and 
of BEC as symmetry breaking pl|j2^ ] . In non-relativistic 
theories these techniques are however nearly unemployed 
and we want to show in the present article that they also 
can be used here in a very effective way. 

The plan of the paper is as follows. First we develop 
the quantum statistics of N nonintcracting bosons in a 
finite cavity of arbitrary shape. We will exemplify the 
use of heat-kernel techniques and of the Mellin-Barnes 
integral representation for the calculation of the parti- 
tion function. We will explicitly show what influence 
the boundary and its shape have on all thermodynami- 
cal quantities, especially on the critical temperature at 
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which BEC occurs. In Sec. Ill we describe an alternative 
approach where sums are converted into ordinary inte- 
grals with an appropriate density of states factor. We 
will show that the density of states is determined by the 
heat-kernel coefficients of the associated Schrodinger op- 
erator. The next section is devoted to the discussion of 
some specific examples like the rectangular box and the 
spherical cavity with various boundary conditions. In 
the final section we present a short summary of our main 
results. 



II. QUANTUM STATISTICS OF A FREE BOSE 
GAS IN A D-DIMENSIONAL FINITE CAVITY 

Let us consider a system of N noninteracting bosons 
in a finite cavity M with boundary dM and with single 
particle state energies E^ determined by the Schrodinger 
equation 



2m 



A<j> N (x) = E N 4> N (x); <t>N\ x< z d M = •••> 



where we do not need to fix here the boundary conditions 
imposed on the field because our treatment will be quite 
general. 

In the grand canonical approach, the partition sum 
then reads 



<l : 



EM* 



(2.2) 



with the fugacity z = exp(/3//), fi being the chemical 
potential and j3 = l/(kT). In the discussion of BEC 
the ground state always plays a special role and for that 
reason we write 



q = 9o 
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(2.3) 
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Here the prime indicates that the ground state is to be 
omitted and q is the contribution of the ground state 
with energy Eq and degeneracy do, 



(2.4) 



For the calculation of the partition sum, equation ( |2.3| 
we will first expand the logarithm to obtain 
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For the evaluation of this kind of expressions it is very 
effective to make use of the Mellin-Barnes type integral 
representations, 



^- J daT(a)v- a , 



(2.6) 



valid for 5Rw > and c € IR , c > 0. Equation (2.6) is 
easily proven by closing the contour to the left obtain- 
ing im med iatel y th e power series expansion of exp(— v). 
Using (2.6) in (2.5) we find (defining \i c = Eq) 
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xV^r / daT(a)((3n)- a (E N -E )- a . (2.7) 



At this stage we would like to interchange the summa- 
tions over N and n and the integration in order to arrive 
at an expression containing the zeta function associated 
with the Schrodinger equation (2.1), 



C( s ) = ]T' (e n - n c y 



(2.8) 
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General zeta function theory tells us that the rightmost 
pole of C(s) is located at s = d/2, see for example [ fj3[ , 
followed by poles at s = (d - l)/2, . . . , 1/2; -(21 + l)/2, 
i € No- In general the rightmost pole appears at s = d/m 
where d is the dimension of space and m is the order of 
the elliptic differential operator, here two. Furthermore 
we need the polylogarithm, 



Li n (x) =E^T' 



i=i 



(2.9) 



basic properties of which may be found in Q,g5|. As 
we will see in the following, the treatment of arbitrary 
dimension d creates no additional complications. But 
it may not only be of academic interest. In this con- 
text we mention the analogies between bosons in a two- 
dimensional box and bosons in a one-dimensional har- 
monic oscillator potential analysed in Ref. [p6| . Similar 
analogies between higher dimensional cavities and exter- 
nal potentials in dimensions d — 1,2,3 remain to be ex- 
plored. 

Due to the above remarks, in order that the summation 
and integration might be interchanged we have to impose 
that 5Rc > d/2 to obtain 
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(2.10) 



This is a very suitable starting point for the analysis of 
certain properties of the partition function q. Closing the 
contour to the right corresponds to the large-/3 expansion; 
closing it to the left to the small-/? expansion. (The rel- 
evant dimensionless expansion parameter will be made 
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cle ar lat er). To the right of the contour the integrand 
in fl2.1C|) has no poles, which means that the large-/? be- 
haviour contains no inverse power in (3. One might show 
however, that the contribution from the contour itself is 
not vanishing at infinity leading to exponentially damped 
contributions for (3 — > oo, the well known behaviour of 
partition sums at low temperature. 

However, the relevant range for the analysis of BEC is 
the small-/? range. (As we will see we are going to obtain 
a reliable expansion if the de Broglie wavelength is small 
compared with the typical extensions of the cavity. For 
/^m-sized cavities this will be true near the transition 
temperature.) Thus we close the contour to the left and 
pick up only the leading two terms at a — d/2 and (for 
d> 1) a = (d- l)/2. We find 

q = q + T{d/2)p- d,2 Li {d+2)/2 (e-^-^) Res C(d/2) 
+r p- {d - 1)/2 L l[d+1)/2 (e-«f-*) 



x Res C((rf- l)/2) 



(2.11) 



with the residues of £(s) denoted by Res £. The case 
d = 1 is special in that ((s) has no pole at s = (d — 
l)/2 = 0. For that dimension the subleading contribution 
contains £'(0) of which, as a rule, no detailed information 
is available. So later, for d = 1, we will consider only the 
leading contribution coming from Res £(1/2) and which 
is correctly given above. 

In order to obtain a reliable approximation (in any 
dimension) at the critical temperature where BEC occurs 
we are thus left with the task of determining the leading 
poles of the zcta function associated with the Schrodinger 
equation. Concerning this determination there are deep 
connections between the zeta function of an operator and 
its heat-kernel defined as [B7] 



K{t) = £' e 



-t(-Ejv-Mc) 



(2.12) 
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The small- i behaviour of (2.12) is of special relevance 
here, 



K(t) 
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(2.13) 
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where for a generally curved spacetime in the meantime 
the first 6 coefficients are known for an arbitrarily shaped 
smooth cavity ]23] and for an arbitrary second order ellip- 
tic operator with leading symbol the metric of the space- 
time. We will need here only the first two coefficients 
which for the Schrodinger operator at hand read 
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(2.14) 
(2.15) 



with V the volume of the cavity, dV the area of its bound- 
ary, and b being a parameter depending on the boundary 
conditions, 

^ ( — 1 for Dirichlet boundary conditions 

1 1 for Neumann boundary conditions. ^ ' ' 

Although to the authors knowledge Neumann boundary 
conditions do not occur in physical applications of non- 
relativistic theories considered here, we include them for 
mathematical completeness and because they do not lead 
to any extra difficulty in calculation. The relevance of 
Neumann boundary conditions in physics stems from the 
treatment of the electromagnetic field in the presence of 
ideal conducting plates, or the bag model in QCD. 



As is easily seen, the residues of £(s), equation ( 2.8| ) 
might be ex press ed through the heat-kernel coefficients 
in equation (2.12). Using the integral representation for 
the T-function we write 



COO 
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dt i«-l e -(^-/*c) < 



(2.17) 



Then we split the integral into t S [0, 1] and t S [1, oo). 
For t he first interval we use the asymptotic expansion 
( |2.13 ), the contribution from the second interval is an- 
alytic in s since the pole terms come from the t — * 
behaviour of the integrand. Thus we find 
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+ finite pieces 
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and read off 

Res C(d/2) = 

Res C((<f- l)/2) 
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(2-18) 

(2.19) 
(2.20) 



Additional poles are located to the left of the above ones. 
The associated residues will depend on the details of 
the boundary of the cavity like its extrinsic curvature 
form. Having however nowadays experimental situations 
in mind, these are probably negligible and so we are not 
going to present more than the subleading order in our 
results. In doing this, results will be very comprehensi- 
ble, but let us stress that higher orders could be included 
easily. 

Usi ng equations ( |2.19| , |2.20|]2~l"4| , |2.15| ) in the result 
( |2.11 ) for the partition sum, we find to the two leading 
orders (as mentioned, for d = 1, now and in the following, 
only the leading piece is correct and will be considered) 
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q = q + Li (d+2)/2 (e-f"' ~' n \ — 
b dV 



V 



with the de Broglie wavelength 

h 



A 7 
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(2.21) 



(2.22) 



Here it is clearly seen, that the expansion parameter is 
given by (typical extension of the cavity)/ At- Our results 
are accurate for large values of the expansion parameter; 
this is, the expansion obtained is an expansion about the 
thermodynamic limit, which is the relevant one for con- 
sidering BEC. It is useful to introduce the dimensionless 
expansion parameter 



(2.23) 



with the definitions V = L d and dV = «£ d_1 . The 
dimensionless parameter k contains information related 
to the shape of the boundary. Then eq. (2.21) reads 



nb . 



-Li 



(d+i)/2 I e 



r 



(2.24) 



At sufficiently low temperatures, where the ground state 
population is large, the approximation fx ~ \x c = Eg 
holds. Then using Li n (l) = ( R (n), n > 1, eq. ( 2.21 ) 
simplifies to 



q = q Q + ( R ({d + 2)/2) 



V_ 



b dV 
- A CR((d + l)/2)^r + 



(2.25) 



which is a very good approximation at some range below 
the critical temperature. 

Let us now continue with the particle number, 



Here one has 

N 



i e 



-0(E N —fi c ) ' 



with the ground state occupation number 
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N = 



,J3(ji c -h) _ i ' 



(2.26) 



(2.27) 



(2.28) 



(Recall that fi c = Eq.) The high temperature expansion 
of N is easily obtained by noting 



= -Li^x). (2.29) 
ox x 

It reads 

N = N + Li d/2 (e-^-^)^d 

+ ^L l{d _ l)/2 (e-^-»))e- 1 + ... (2.30) 

The critical temperature T c (in the absence of a phase 
transition) can be defined as the temperature at which 
the number of particles in the ground state begins to be- 
come large. We will set N — fNo with / representing 
the fraction of particles in the ground state. If we as- 
sume iVo >> 1, then for temp eratures close to T c we 
expect /i ~ /i c . (S ee E q. ( 2.28 ).) The behaviour of the 
polylogarithms in ( |2.30| ) d epen ds on the dimension d. For 
d > 3 all of the terms in ( 2.30| ) are finite as /i — > [i c , and 
we can approximate ( |2.3C ) by 



N ~ { R (d/2)^ 



(2.31) 



The critical temperature is defined via £ c obtained by 
taking T = T c in ( |2.23[ ), This assumes / << 1. Defining 
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N 
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(2.32) 



which gives the critical temperature in the bulk (or ther- 
modynamic) limit, 



Tn = 



N 



2irmk \VC,R{d/2) 



2/d 



(2.33) 



and assuming that £ c = Co(l + 7), with 7 -C 1, to leading 
order we can approximate (we use the temperature here) 



r, = r 1- 



Kb ( R ((d -l)/2) 1 
2d c, R d - 1)/d (d/2) N^ d 



(2.34) 



It is seen, that depending on the boundary conditions im- 
posed, the critical temperature can increase or decrease. 
The magnitude of the correction behaves like _/V _1 / d and 
is typically not negligible. It is seen explicitly that the 
corrections are going to vanish for N — > 00. 

Let us now consider the lower dimensional cases and 
we start with d — 3. Then we need 

Ln (V<3(^-^ = -ln(l - (e^^)). (2.35) 



Near /x ~ /i c we use eq. (2.28) to approximate 

1 \ 1 



p(ji c - n) = In 1 + 



N 



N ' 



(2.36) 



Then in 3 dimensions the analogue of eq. (2.31) reads 
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N(l - f) 



(2.37) 



Here it is not possible to put / = in order to define a 



critical temperature, since in orde r to o btain (2.37) we 
have assumed Nq >> 1. (See Eq. (2.36).) This roughly 



reflects the fact that in 3 dimensions (as well as in 1 and 2 
as seen below) the number of excited states is lower than 
in higher dimensions, and as a result, in the temperature 
range considered, particles have to reside in the ground 
state in order that the thermodynamic equation for N 
is fulfilled. However given / ^ and N, this might be 
solved for the temperature, 



T = T {l--(f + a) 



with 



and 



To 



N 



■2tt ink Vu,T!/2)y) 



nbln(fN) 

4C^ /3 (3/2)JVV3' 



2/3 



(2.38) 
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valid near fi ~ /x c . This parallels completely eq. (2.34). 
Defining the critical temperature T c as the temperature 
where a given fraction / of the number of particles N 
is in the ground state, eq. ( 2.38 ) can easily be used to 
determine T c . This equation also shows that in the ther- 
modynamic limit, / may be put to zero and To is a "real" 
critical temperature. This reflects the fact that in three- 
dimensional unbounded space the ideal Bose gas exhibits 
Bose-Einstein condensation in the sense of a phase tran- 
sition. 

We proceed with d — 2 and need in addition 



Li 1/2 (e 



(2.41) 



valid for /i ~ /i c . Using the same approximation as be- 
fore, (3(nc — m) — V-^b, this results in 

N(l-f)=eHfN) + ^^^jN+... (2.42) 
Solving to leading order for £ given / ~ we find 



or in the temperature 



T, 



N 



HfN) 
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Finally in d — 1 considering only the leading contribution 
results in 



with 



T c = 



h 2 



N 
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(2.45) 
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In dimensions d — 1 and d = 2 it is clearly seen that 
we cannot put / — 0. This reflects the fact that even in 
the thermodynamic limit Bose-Einstein condensation (as 
a phase transition) is not possible. 

Another interesting quantity of the system is its inter- 
nal energy defined by 



U = 



d fi d 



(2.47) 



Of course the energy of the system is the sum of the 
single particle energies En multiplied by the occupation 
number, 



U = U 



N 



-/3(B N -fla) 



Here Uq is the zero-point energy, 
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U = 
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To the approximation considered, U is most easily found 
by using (2.11) and ( p. 29 ). As an immediate result we 
have 
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dq 
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With the already presented expansions around [i ~ \i c 
the behaviour near the phase transition is easily dis- 
played. 

Finally we are left with considering the specific heat, 



C 



dU 
df 



(2.51) 
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The slightly more difficult point here is, that N has to 
be considered as fixed and so (x has to be considered as 
a function of iV and j3. We will need (djd fi) (f3fi), which 
is found by differentiating equation (2.27) given that N 
is fixed, 
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dN 
~d[3 



= 0. 



(2.52) 
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Using eq. ( 2.3C ) to the leading order, condition ( 2.52 ) 
yields 



|[/^o-M)] 



-/3(Hc — A 1 ) 



(2.53) 



It is seen, that in the approximation // = fi c one finds 
the above equation to yi eld z ero as it should. Neglecting 
pieces coming from eq. ( [2.53 ) (one can show that this is 
justified for d < 7) we only have to use ( [2. 50 ) to find 
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(2.54) 



showing once more the decrease or increase in the spe- 
cific heat depending on the boundary conditions. Using 
previous formulas, expansions for /1 ~ /i c and for specific 
dimensions are easily found. 

A slightly different treatment, which we do not pursue 
in this paper, could also be used. This uses the effective 
fugacity (for the anisotropic harmonic oscillator see |14|) 



(2.55) 



where E\ is the first excited level with degeneracy d\. 
Whereas in the calculations presented above only the 
ground state has been treated separately. We can now 
separate the ground state and the first excited level to 
find 

V 

q = qo+ diLitiZeff) + Li {d +2)/2(z e ff)^j 

Arp 

b dV _. 

+ lJdrT L Hd+l)/2(Zeff) +■■■ 

This expansion is a very good approximation even below 
the critical temperature (for the harmonic oscillator see 
PI). All thermodynamical properties are obtained in 
the same manner as before. 



III. DENSITY OF STATES METHOD 

An alternative approach to the use of the Mellin- 
Barnes contour integral representation is to try to convert 
the sums for the thermodynamic quantities directly into 
ordinary integrals with an appropriate density of states 
factor. This has been widely used in a variety of problems 



in statistical mechanics. ( See Re f [p9[ for example.) Re- 
cently Grossmann and Holthaus [|10|Jll| have used such a 
density of states method to study the statistical mechan- 
ics of particles in a harmonic oscillator potential trap, and 
in a 3-dimensional cubic box fil| with Dirichlet boundary 
conditions. The harmonic oscillator trap is characterized 
by a density of states 



E 2 E 



(3.1) 



where Q = {uji^loz) 1 ^ is the geometric mean of the 
anisotropic harmonic oscillator frequencies. In Ref. ]3C| ] 
we showed how it is possible to evaluat e 7 analytically, 
and obtain additional corrections to (pTl|). (See also 
fl4j|.) The purpose of the present section is to discuss 
the method in greater detail in a more general setting, 
and to show how results equivalent to those of the pre- 
vious section may be obtained with the density of states 
method. The basic idea behind the method can be found 
in Ref. ]3l| ] with a number of examples illustrated. 

Suppose that we have a self-adjoint differential opera- 
tor A which possesses a set of eigenvalues {A at} which are 
all positive. Here N represents a multi-index which labels 
the eigenvalues. For the case of statistical mechanics in 
a cavity, A can be taken to be the Hamiltonian with Xn 
the energy levels Em- The eigenvalues depend on the de- 
tails of the cavity and boundary conditions imposed, but 
the method we describe can be used regardless of such 
details. Let A/"(A) be the number of modes for the bound- 
ary value problem in the cavity for which the eigenvalues 
A at < A. We can write 



N{\) =^0(A-Aat) 



(3.2) 



where 9{x) is the Heaviside distribution (or step function) 
defined by 



(3.3) 
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= 



6(x) = 
This gives the result 



(3.4) 



A N <A 



Ajy— A 



as used by Baltes and Hilf . 

The aim now is to treat the eigenvalues as a continuous 
distribution rather than a discrete set by introducing the 
eigenvalue (or spectral) density p(A). We define 



p(X)d\=JV(X + d\) -JV(X) 
dM(X) 



dX 



-dX 



(3.5) 
(3.6) 



Use of ( [3.2] ) for A/"(A), along with the distributional iden- 
tity 6'(x) = 5{x) gives 
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P(A) 
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S(X-X N ) 



(3.7) 



If we use the familiar exponential representation for the 
Dirac delta distribution, then 



P(A) 
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(3.8) 



We may now translate the integration contour in (3.8) to 
the right by an amount c where c G fft with c > 0, and 
assume that we are permitted to interchange the order 
of the summation and integration. This results in 



P(A) = 
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dte tx K(t) , 



where 



tA? 



(3.9) 



(3.10) 



Eq. (|3.9[ ) shows that p (A) i s the inverse Laplace transform 
of K(t). Inversion of (3.9) results in 



K(t) 



d\e- tx p(X) 



(3.11) 



Another way of arriving at ( p.9h is t o note that ( 3.11 ) 
follows directly from ( |3.7| ) and ( |3.10| ), establishing that 
p(X) is the Laplace transform of Kit). The inversion 
formula for Laplace transforms gives ( |3.9[ ) immediately 
without the intermediate steps in ( |3.8j ). 

The po int o f (3.9) and (3.11) is that given K(t) as de- 
fined in ( 3.10 ) we can evaluate the density of states in 
a straightforward way. In certain simple cases, such as 
the harmonic oscillator potential, where the eigenvalues 
can be writ ten d own explicitly, it is possible to perform 
the sum in (3.10) and evaluate K{t) in closed form. (See 
Ref. ||(J for a case of this.) In the generic case, where 
the eigenvalues are not known explicitly, we can still ob- 
tain a result for the density of states by using the known 
asymptotic expansion for Kit) as t — > p7| . 

Assume that as t — > we have the generic expansion 



K(t) ~^2 c ^ r ' + °( t ~ rk+1 ) 



(3.12) 



for some coefficients Cj and powers with n > > 
•__> rfc. The simplest way to evaluate p(X) is to use 
( |3.11 ) and the relation 



1 



dX A 2 ~ 1 e~* A 



r(*) Jo 

which is valid for 5ftz > 0. It is now easy to see that 



(3.13) 



p(x) - E 



1 nn) 

Integrating this result with respect to A results in 



(3.14) 



Af(X) ~ 



c, 



~i r(r< + 1) 



(3.15) 



Although the analysis we have presented makes no at- 
tempt at proper mathematical rigour, the main result in 
Q3.14 ) of ( |3. 15 ) is in agreement with a refinement of Kara- 
mata's theorem due to Brownell |32|. (A nice account is 
contained in Ref. It is important to emphasize that 

the result for p(A) or Af(X) refers really to the average 
number of eigenvalues. Generally the number will fluc- 
tuate around some average value. These fluctuations can 

Another point 



51 



33 



be important in some contexts 
wo rth st ressing is that our results have assume d tha t 
in ( 3.12| ) is positive. The results ( |3.14| ) and ( 3.15 ) are 
actually true for all if we define l/r(rj) = when 
Ti = 0,-1,-2,.... The establishment of this result re- 
quires a more powerful approach than Laplace transforms 



The density of states (3.14) is determined by a knowl- 
edge of the coefficients c; appearing in the asymptotic 
expansion of K(t). These coefficients can be evaluated 
from a knowledge of the generalized ^-function associated 
with the eigenvalues Ajv as discussed in the previous sec- 
tion. Defining 



c(*) = E A 



— s 
N 



(3.16) 



the analysis in Sec. 2 is easily modified to show that £(s) 
has a simple pole at s = ri with residue 



ResC(s)| s 



r(r0 



for n 7^ 0,-1, —2, . . .; and 

C(s = -r i ) = (-ir(-r i )!c i 



(3.17) 



(3.18) 



for n = 0,-1,-2,.... Therefore the coefficients ap- 
pearing in the density of states may also be found from 
a knowledge of the generalized ^-function. This provides 
a direct link between the density of states approach and 
the method described in Sec. 2. 

We will conclude this section by showing how the den- 
sity of states method may be used in the specific case of 
an ideal gas confined in a general cavity. The g-potential 
defined earlier was written as 



q = Qo + Qe 



where 



qo = —do In ^1 — ze 



PE 



(3.19) 



(3.20) 



7 



was the ground state contribution, and 



For periodic boundary conditions, 



22, e nP(p-ii c ) 



E 



E' e ""' 



0{E N -E o ) 



(3.21) 



represents the contribution from the excited states. The 
sum over energy levels in ( pi] ) can now be converted 
into an integral over a continuous energy vari able b y in- 
troducing a density of states factor. From ( 3.14 ) and 
using the coefficients Cj found in Sec. 2 we have 



p{E) 



A£r(d/2) + 1 x^rdd- 1)/2) 



(3.22) 



if we keep only the first order correction to the bulk ex- 
pression as before. To this order we may approximate 



^ e n/3{p-fi B ) 



fe ~V- / dEp(E)e 

n=l U J Ei-Bo 



-T10E 



(3.23) 



Replacing the lower limit in the integration with 0, we 
obtain 



b ^-(d-l 



\dV)Li {d+1)/2 {e-^-^) (3.24) 



in agreement with our earlier method. The replace- 
ment of the lower limit on the integral with assumes 
f3E << 1 so that the temperature is well above the 
temperature associated with the ground state energy (i.e. 
kT >> Eq.) If the energy gap Ex—Eq is of the same order 
of magnitude as the ground state energy, then {3E± < < 1 
as well. Results for the particle number and other ther- 
modynamic quantities follow as before. 



IV. APPLICATION OF THE RESULTS TO SOME 
CAVITIES 

This section will discuss some specific examples like 
the rectangular box and the spherical cavity with various 
boundary conditions. The rectangular box has been the 
subject of considerable investigation. (See Refs. [|||J35) 
for reviews.) In order to get all thermodynamic quanti- 
ties as derived in Sees. II and III, we only need to give the 
constant &, the volume of the cavity V, and the bound- 
ary dV. Another interesting quantity is the lowest energy 
eigenvalue as the critical value for the chemical potential. 



A. Rectangular box 

For the rectangular box with side lengths L\, . . . , L d 
one can consider for example periodic, Dirichlet and 
Neumann boundary conditions. In cartesian coordinates 
x = (x\, . . . ,Xd) with Xi e [0,Lj]. 



<j>(xi, ...,x l ,...,x d ) = 4>{xi, ...,Xi+Li,.. .,x d ), (4.1) 
with i = 1, . . . , d, the energy eigenvalues are 

2 



E ni ,...,n d — E ( r . 



2?r 



rii G Z. 



(4.2) 



Periodic identification means, that effectively there is no 
boundary. (Equivalently we can regard V as the e?-torus.) 
In this case b — 0. Furthermore, one has for the volume 



V d = f[L i 



(4.3) 



and [i c = Eq = 0. 

For Dirichlet boundary conditions, 

4>(xi, . . . ,Xi_i,0,Xi+i, ...,Xd) = 0, 
4>(xi, . . . ,Xi-i,Li,Xi+t, . . . ,Xd) = 0, 

one gets 



d / \ 2 
nut 



(4.4) 



As mentioned, see eq. ( p.lq ), b = — 1, the volume is as 
before and furthermore 



07 = 2£vj« 



(4.5) 



i=i 



(i) 

Here, V d _ 1 is the volume of the i'th side of the rectangular 



cavity, V d } l = L\ . . . Li—xLi+x . . . L d ■ Furthermore, 

d / 2 

L~ 



(4.6) 



If we specialize to the d-cube, then dV = 2dL d ~ 1 which 
gives k = 2d. 

Finally, for Neumann boundary conditions, 

(d/dxi)<f>(xi,. . . ,Xd)\ Xi =o = 
(d/dxi)(f>(xi, . . .,x d )\ Xi = Li = 0, 



the energy eigenvalues read 



E n 



E 



rii e N , (4.7) 



with /u c = Eq = 0, b = 1 and V, dV, as before. 

All thermodynamical quantities may now be given im- 
mediately with the formulas in Sec. II. Especially for the 
critical temperature one has that periodic boundary con- 
ditions give a critical temperature identical to the ther- 
modynamic limit (a result of being no boundary there), 
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Dirichlet boundary conditions increase it and Neumann 
boundary conditions decrease it. 

Looking at the distribution of the eigenvalues around 
the ground state of the system an intuitive explanation 
of this behaviour is possible. Let us assume for this 
equal compactification sides Li. It is clear, the lower the 
density of eigenvalues near the ground state, the higher 
the transition temperature will be. Then the following 
might be stated. Neumann boundary conditions lead to 
the lowest condensation temperature because it has the 
smallest spacing between the ground state and the first 
excited level. This compensates the doubled number of 
eigenstates in the periodic case, the eigenvalues which 
are however four times higher. The difference between 
Neumann and Dirichlet boundary conditions is based on 
the same observation; in addition the degeneracy of the 
eigenvalues is the same for this two boundary conditions. 
The difference between periodic and Dirichlet boundary 
conditions is not obvious, however the results show, that 
the slightly smaller level spacing for Dirichlet boundary 
conditions is more than compensated by the doubling of 
the degeneracy for periodic boundary conditions. 



B. Spherical cavities 

For spherical cavities an explicit knowledge of the 
eigenvalues is not at hand and for that reason also a nu- 
merical treatment of the partition sum is not immediate. 
For this kind of examples our analytical approach is most 
useful. 

It is convenient to introduce a spherical coordinate ba- 
sis, with r — \x\ and d — 1 angles f2 = (6*i, . . . , 8d-2, <p)- 
In these coordinates, a complete set of solutions of the 
Schrodinger equation together with one of the mentioned 
boundary conditions may be given in the form 



n d/2 R d 2n d / 2 R d - 1 
V = ov 

r(d/2 + i) ' 



(4.10) 



(f>Lm,n(r, CI) = r 1 2 J 



^Wi, n r/h)Y l+sL {tt), (4.8) 



with J; + ( £ ;_2)/2 being Bessel functions and lj+d/2 hy- 
perspherical harmonics |36|. Here, Ei >n = wf n and the 
wi,n(> 0) are determined through the boundary condi- 
tions by 



^2mwi, n R/H) = 0, 
for Dirichlet boundary conditions and 



l-d/2 , , — 

= — J l ! d-2 (y 2mwi.n.R/ n) 



(4.9) 




J' d -2 (V2mwi, n r/h) \ r =R 



for Neumann boundary conditions. Although none of 
these eigenvalues is known explicitly, in order to obtain 
an approximation of all thermodynamic quantities, it is 
already enough to state the volume and the boundary 
volume of the spherical cavity, 



T(d/2) ' 

and to apply the results of the previous sections. 

V. DISCUSSION AND CONCLUSIONS 

We have described two methods for treating the ideal 
Bose gas using the grand canonical ensemble in a general 
cavity. The first method, described in Sec. II, treats the 
partition function as a sum over the discrete energy levels 
by making use of the Mellin-Barnes integral representa- 
tion for the exponential function. The second method 
shows how to obtain an adequate density of states fac- 
tor so that the sum over discrete energy levels can be 
replaced with an integral. Both methods make use of 
the so-called heat-kernel coefficients for the Schrodinger 
operator, and equivalence between the two methods is 
shown in quite a general setting. 

Our analysis shows clearly the way in which the ge- 
ometry of the cavity enters in all thermodynamic prop- 
erties. In the limit considered, namely £ = L/Xt 1, 
the leading correction to the bulk limit term £ d enters 
via the property k — (dV)/V( d ~ x " d . Thus the leading 
order correction is completely described by the volume 
and area of the cavity, with all finer details washed out 
for £ 3> 1. Finer detail would be present in the next 
to subleading order in terms of the extrinsic curvature 
of the boundary. (The extrinsic curvature describes how 
the normal vector to the boundary varies when moving 
along the boundary.) 

These ideas are as well applicable to the canonical en- 
semble. Furthermore, due to the connections between the 
microcanonical and the grand canonical approach stated 
recently in |]37|"|39|| it seems possible to develop this ap- 
proach also for the microcanonical ensemble. 
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